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ABSTRACT 

The effects that farge scale fluctuations had on small scale isothermal modes at the 
epoch of recombination are analysed. We find that: (a) Albeit the fact that primordial 
fluctuations were at this epoch still well in the linear regime, a significant nonlinear ra- 
diation hydrodynamic interaction could have taken place, (b) Short wavelength isother- 
mal fluctuations are unstable. Their growth rate is exponential in the amplitude of the 
large scale fluctuations and is therefore very sensitive to the initial conditions, (c) The 
observed CMBR fluctuations are of order the limit above which the effect should be 
significant. Thus, according to their exact value, the effect may be negligible or lead to 
structure formation out of isothermal fluctuations within the period of recombination, 
(d) If the cosmological parameters are within the prescribed regime, the effect should 
be detectable through induced deviations in the Planck spectrum, (e) The sensitivity 
of the effect to the initial conditions provides a tool to set limits on various cosmologi- 
cal parameters with emphasis on the type and amplitude of the primordial fluctuation 
spectrum, (f) Under proper conditions, the effect may be responsible for the formation 
of sub-globular cluster sized objects at particularly high red shifts, (g) Under certain 
circumstances, it can also affect horizon sized large scale structure. 
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1 INTRODUCTION 

The problem of structure formation in the Universe has 
probably been one of the foremost and most studied question 
in cosmology. Perhaps the greatest achievement of cosmol- 
ogy was the prediction and following discovery of the Cos- 
mic Microwave Background Radiation - the CMBR (Pen- 
zias & Wilson, 1965) and its anisotropy (Smoot et al. 1992). 
Not only did it provide support for the Big Bang theory, 
it proved that structure is a result of small fluctuations 
growing into large inhomogeneities and not vice versa. The 
fact that the observed fluctuations in the CMBR are small 
(AT/T ~ 10 -5 ), naively implies that one can treat the fluc- 
tuation within the linear approximation when modeling the 
evolution of structure before the time of last photon scatter- 
ing, i.e., one can assume that fluctuation modes of different 
wavelengths are completely decoupled from each other be- 
fore radiation-matter decoupling. 

The first to develop a linear theory for the perturbed 
Fridmann - Robertson - Walker metric was Lifshitz (1946), 
and it can be found in various text books such as Weinberg 
(1972), Peebles (1980) and Kolb & Turner (1990). The the- 
ory has since then been applied to describe cosmologies with 
various components and various initial parameters. Baryonic 
matter, radiation and other massless particles, cold (mas- 



sive) or hot (light) dark matter are some of the ingredients 
that enter the primordial soup. While other parameters such 
as the Hubble constant Ho, the vacuum energy through the 
cosmological constant A, and the initial spectrum affect the 
evolution. A review of various current cosmological models, 
their evolution and implication can for example be found in 
White et al. (1994) with emphasis on the microwave back- 
ground radiation and in Primack (1997). Less current re- 
views that cover the basic principles of structure formation 
are found in the aforementioned textbooks. 

One of the major parameters that affect the qualita- 
tive behaviour of our universe and is highly relevant to the 
present paper, is the type and form of the primordial spec- 
trum of fluctuations. The fluctuations are generally classified 
into curvature (or adiabatic) and isocurvature (or isother- 
mal) fluctuations. The first arise naturally in inflationary 
scenarios (e.g., Liddle and Lyth, 1993, and ref. therein). 
They are fluctuations of both the baryonic fluid and the ra- 
diation and they propagate at the adiabatic speed of sound 
that is equal to c/v3 if the radiation energy density domi- 
nates. These waves are found to decay on scales smaller than 
the Silk scale, namely, scales comparable to galaxy sized ob- 
jects (Silk, 1967). Therefore, a top-down structure formation 
is a natural consequence of adiabatic perturbations. Small 
scale objects can form after recombination (and initiate a 
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bottom-up scenario) only if one adds the undamped per- 
turbation of a cold component as is the case in the stan- 
dard CDM model (Peebles 1982, Blumenthal et al. 1984 
and Davis et al. 1985) and its variants. Isocurvature (or 
isothermal) fluctuations on the other hand, are a natural 
consequence of topological defects formed in phase transi- 
tions (e.g, strings, monopoles and textures) or if more than 
one field contributes significantly to the energy density dur- 
ing inflation. They correspond to fluctuations that alter the 
entropy density but not the energy density. Unlike the first 
type of fluctuations, these do not suffer from Silk damp- 
ing. Consequently, objects with a mass as small the post- 
recombination Jeans scale, namely the size of globular clus- 
ters, can form after decoupling. 

It is interesting and important to note that even if the 
primordial spectrum was purely adiabatic, isothermal per- 
turbations of a given wavelength are formed as the second 
order effect of Purcell clustering from adiabatic waves of 
similar wavelengths (Press & Vishniac, 1979). Through this 
effect, nonintcracting particles (baryons) that are viscously 
coupled to a stochastically oscillating background gas (ra- 
diation fluid) undergo secular clustering. Originally, it was 
thought that the effect can produce a bottom-up scenario 
even from a pure baryon model with an adiabatic spec- 
trum. However, the typical post-recombination Jeans scale 
amplitude of the isothermal waves will only be pi SO /p ~ 
(1CP 2 — 1) x (Sp ad /p) 2 , with Sp ad /p the typical adiabatic 
fluctuations at horizon crossing. Namely, if the primordial 
spectrum is flat and adiabatic, the typical isothermal fluc- 
tuation at recombination is roug hly 8p lso /p~ 1(T 10 -Hn 8 . 
For a flat isothermal primordial spectrum, one should expect 
typical amplitudes of pi SO /p ~ 10~ 5 — 10~ 4 . 

The smallness of nonlinear effects such as the Purcell 
clustering (or shock waves which are of an even higher order) 
led to the consensus that the evolution of the fluctuations 
can be treated linearly and modes of different wavelengths 
are decoupled from each other. This delays the nonlinear 
treatment to the time when radiation does not play any 
dynamic role anymore. At face value, it certainly appears 
to be the case as 5p/p, v/c, and ST/T are all much smaller 
than unity. One should nevertheless be extremely careful 
when assuming linearity, especially in view of the fact that 
not all of the dimensionless parameters are actually smaller 
than unity. 

One of the dimensionless numbers that appears in the 
solution of the nonlinear fluctuations' equations of motion 
and that is not small at all is the root of the radiation to 
ga s pressure ratio. Just before recombination one finds that 
\/ Prad/ Pgas ~ 10 5 ! Consequently, we expect that the inter- 
action between the radiation and matter at this period will 
have profound impact on the evolution of the fluctuations. 

In this paper we examine the linear hypothesis and its 
validity by adding the force large scale perturbations exert 
on short wavelength waves. We begin in §2 by overviewing 
the problem of solving the nonlinear equations of motion 
and estimating the effect with a very simple analysis. We 
proceed in §3 to write the Newtonian equations describing 
the evolution of short wavelength isothermal modes. In §4 
we analyse the simplified solution, while in §5 we proceed 
to estimate the effect in a few cosmological scenarios. In 
§6 & §7, we study the possible ramifications to structure 
formation and study the possibility of measuring and using 



this effect in the study of cosmological parameters. In §8 we 
show that the effect can influence large scale structure as 
well. 



2 THE EFFECT 

At the lowest order of approximation, waves of different 
wavelengths are decoupled from each other. This suggests 
that a small scale isothermal wave will, at this order, witness 
an isotropic and homogeneous medium around it. However, 
at the next order of approximation, one has to include the 
small perturbations (of order Sp/p ~ 1CP 4 ) of all the other 
scales in the vicinity of the isothermal wave. 

One can simplify the problem by assuming that the 
isothermal wave is localised to a region of size d that is 
at least several times the wave's wavelength (such that the 
wave parameters are well defined). The scale d separates be- 
tween two types of perturbations: those with a wavelength 
A < d and those with A > d. The contribution for example 
to the radiative force of the former group vanishes in the 
lowest order as the spatial correlation between the wave and 
perturbation yields a net vanishing result. The latter group 
can contribute a net effect only if the temporal correlation 
between the perturbation and the isothermal wave does not 
vanish as well; namely, it will do so only if the perturba- 
tion does not have enough time to oscillate rapidly inside 
the region under consideration. This implies that for a given 
period At, the contribution of waves with uAt i> 1 aver- 
ages out leaving no net contribution. The only waves that 
have a net contribution are those for which to At 1; that 
is, if d is chosen smaller than the wavelengths of waves with 
lo ~ 1/At, its exact value is unimportant. Hence, we solve 
for short wavelength isothermal waves in a region where the 
long wavelength perturbations can be considered constant 
and homogeneous. We cannot however, assume that the lo- 
cal environment is isotropic. 

We assume for simplicity that the large scale perturba- 
tions are optically thick adiabatic waves of the form: 

— = 8a sin (k n • x - oj a t) , (1) 
P 

where Sp a is the perturbation of the matter density. It is 
proportional to the temperature perturbations: 

3ST a /T = Sp a /p, (2) 

namely, the fluctuations in the radiation field and in the 
matter are synchronised. The sole force responsible for ac- 
celerating the matter is the radiative force. If observed from 
the baryon rest frame, the force is due to a net radiative 
flux originating from an anisotropy in the radiation field. 
This anisotropy or shear between the radiation and matter 
fluids can be found through the baryon acceleration. 

Using the adiabatic speed of sound «( a ) to relate k a with 
Wo and the continuity equation to relate v a with 8p a , the 
velocity of the baryonic fluid in the adiabatic wave is found 
to be: 

v a = Sv(a) sin (k a ■ x - u) a t) n a , (3) 

with n a being a unit bector in k a 's direction. The acceler- 
ation or the force per unit mass is simply dv a /dt, conse- 
quently, if the opacity per unit mass is k, the net radiative 
flux H is given by: 
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H 



aa 

K 



1 d\g 

k dt 



■ cos (k a • x — uJ a i) n a 



(4) 



Si sin (ki ■ x — u>ii) . 



(•») 



The form for an isothermal wave travelling in a small 
region perturbed by the wave above is: 

5 pi 
P 

with STi = as it is isothermal. It will witness a constant ra- 
diative flux H given by eq. |^ as it is optically thin. However, 
if k is a function of density, the opacity and therefore the 
force per unit mass, become a function of the wave's phase. 
The predominant source of opacity before and during recom- 
bination is Thomson scattering off the free electrons. Con- 
sequently, the only period during which isothermal waves 
vary the opacity is during recombination, when the number 
of free electrons is proportional to p~ 1//2 . Hence, the total 
acceleration is: 

a = kH = k H + SkJI = a a + a; = a a ( 1 — ~ [ , (g) 

V 2 o ) 

The first term results with the large scale acceleration al- 
ready accounted for in the large scale fluctuation. The sec- 
ond term however, incites a force that varies synchronously 
with the small scale isothermal waves. It leads to an insta- 
bility. 

The total power input per unit mass by the latter term 
into the wave is: 

Prad = 3i ■ Vj = \-V {a) S a UJ a COS (ka ■ X - LU a t) fi a 



• Si sin (ki • x — LUit) Viiii 
1 



V( a )LJ a 6 a COS (k a 



x sin (ki ■ x — ujit) 2 v^Sf cos ( 



(7) 



where fii denotes a unit vector in ki's direction and 9 is the 
angle between fii and n a . Again, we related to pi through 
the continuity equation and the isothermal speed of sound 
V(i). When averaging the force over time one has to recall 
that the large scale perturbations are assumed not to vary 
over the relevant period - the duration of recombination. 
Now we see why. Oscillating large scale perturbations will 
average to a zero net contribution. Those that oscillate on a 
longer time scale can be assumed constant: 



5 a cos (k • x — cut) 



(8) 



Using ^cos 2 ^ = 1/2 to average the isothermal oscillations, 
we find: 

1 2 

(jPrad) = -VaLjS^V^Si COS 9 . (9) 

The total energy per unit mass in the acoustic wave is 

e = itf«f /2. (10) 

The e-fold growth rate of the wave's amplitude is therefore: 



r — -^ J - — - — ; „ „ — cos ( 



2e 



2U 2 ' 



1 x V ^ a 

— LUOrt, COS 9. 

4 t) (i) 



(11) 



Let the duration of recombination be At. The growth factor 
G, or the number of e-folds the wave will grow if traveling 
in the direction of the large scale shear (i.e., if 9 = 0), is: 

G = rAt = ^'HiSl ( wA t) = Go-^- with Go = 0(1), (12) 

4 Vtf) 10 



where we have inserted the ratio of the speeds of sound at 
recombination assuming radiation dominance over baryonic 
matter. We have also taken uiAt ~ 1, as waves with a larger 
u> average out. Go is a constant of order unity which cannot 
be found in this approximate estimate. 8$ at recombination 
will have a distribution of values. If however, the typical 
value is of order 10 -4 , then the typical growth rate is Go e- 
folds. Namely, we expect that the nonlinear effect is at least 
as important as the linear evolution of the isothermal wave. 
It is now also evident why isothermal waves are unstable 
while adiabatic are mostly stable. One only has to replace 
Vu) by V( a -) to find that the adiabatic waves are unstable 
only for Ss of order unity and not much less. 

Eq. [13 is useful as an estimate, but a more general ex- 
pression for the growth parameter G should not include an 
assumption on the type, speed, or optical depth of the large 
scale waves, but instead use directly the radiation flux ob- 
served in the baryon rest frame. By using eq. ^| and letting 
kH vary during recombination, we have: 



G 



4 c 



W Jo 



nHdt. 



(13) 



Note that H is the flux observed in the baryon rest frame. 



3 NEWTONIAN EQUATIONS OF MOTION 

The general nonlinear relativistic equations describing the 
evolution of the perturbations are extremely difficult, if not 
impossible, to be solved analytically. We therefore simplify 
the problem until it becomes amenable to an analytical 
treatment without losing the basic physical ingredients. It is 
generally not possible to decouple the system into different 
wavelengths and solve each one separately. Nonetheless, if 
the nonlinearity is small, namely, the interaction between 
modes of different wavelengths is small, we can write the 
equations describing each mode as decoupled equations to 
which a coupling correction is added. Moreover, we have 
seen in §2 that the coupling is essentially an effect large 
modes (adiabatic or isothermal) have on small isothermal 
waves. That is, the large horizon sized modes are largely un- 
affected and can be integrated as uncoupled modes, whereas 
the small modes are mostly affected only by the large modes 
and not by themselves. Consequently, the short wavelengths 
modes can be solved independently as the equation for each 
one is the basic uncoupled linear equation to which a linear 
source term is added. This term depends only on the history 
of the long wavelengths modes. Actually, this statement is 
true only as long as the amount of energy transferred into 
the short wavelengths is negligible when compared to the 
energy in the larger scale. This will be further discussed in 
§8. 

To summarise, the assumptions made in order to sim- 
plify the equations are: 

• The short wavelength fluctuations are pure isothermal 
waves - i.e., there are no perturbations to the radiation field. 

• The fluctuations are of a wavelength much smaller than 
the Horizon size and are therefore solved under the Newto- 
nian approximation. 

• The short wavelength fluctuations are affected only by 
large wavelength modes. 
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Large Scale: 



shear = net radiation flow 




Small Scale 



shear 



(isothermal fluctuations) b 

Figure 1. The physical system we solve. The equations describe 
the evolution of small scale perturbations embedded in an envi- 
ronment that is inhomogeneous on large scales and can have both 
a perturbation in the baryon density S), and in the radiation en- 
ergy density <5 r . The large scale inhomogeneities are the Horizon 
scale perturbations. The small scales are very short wavelength 
isothermal perturbations solved locally in a small region where 
the large scale can be considered uniform, or in other words, that 
the large scale flux is spatially constant. This approximation is 
possible because the small scale wavelengths considered are much 
smaller than the inhomogeneities affecting them. 



with v m being the material velocity. Either way, the total 
radiative force per unit mass f ra d can be written as: 



trad — Hk — 



(16) 



We will work with the latter expression that uses the effec- 
tive velocity of the radiation fluid v r as it is easier to work 
with when coupling the baryons to the radiation. For non- 
relativistic velocities and Thomson scattering, the term can 
be written as: 



c 4 4 v m - v r n e 

trad = —aral 

3 c pt 



(17) 



n e is the number of free electrons, a denotes the radiation 
density constant, and or is the Thomson cross section. 

Normally when linearised and Fourier transformed, this 
term which describes the radiation drag has a contribution 
only from the wave for which the linearized equations of 
motion are written. However, the contribution of waves of all 
scales should be taken into account if through the nonlinear 
interaction they produce a correction synchronised with the 
waves. When linearising, we find the following terms: 



frarf 



(18) 



• The evolution of the large wavelength modes is decou- 
pled and is assumed to have been solved a priori. 

• The short wavelength modes are solved in localised re- 
gions where the effect of the large scales can be considered 
as homogeneous. 

The system solved is described in figure jjj. 

We start by writing the equations of motion describing 
Newtonian isothermal fluctuations. The simplest equations 
describing matter fluctuations of scales much smaller than 
the horizon size and velocities much smaller than the speed 
of light are the Newtonian equations of motion, corrected 
to include the expansion of the universe. These are (see for 
example Kolb & Turner 1990): 



o a a. . „ 
pi + 3-pi + -(r ■ V)pi + p V ■ vi 
a a 

a a _ , „ 

Vl + -Vl + -(r ■ V)Vl -| Vpi + V01 - fracl 

a a po 



4nGpi 



(14) 



where the various quantities have their usual meaning. Here 
v a is the isothermal speed of sound. The last term in the 
second equation is the force per unit mass acting on the 
baryonic fluid. It can be expressed with the help of the flux 
H observed in the baryonic rest frame and the opacity per 
unit mass k. It can also be expressed with the radiation 
"velocity" v r - the velocity with which the baryons need to 
move in order for them not to observe a net radiation flow, 
and a characteristic time r e7 for the transfer of momentum 
between the electron-proton fluid and the photons. If the an- 
gular dependence of the photon energy density distribution 
Sp r (9,x) is expanded into Legendre polynomials: 



Sp r (9,x,t) 

pr(x,t) 



^2(2£ + 1) Pe{cos9)a e (x,t), 



(15) 



then the radiation "velocity" is given by v r = ci/4 and 
the flux observed in the baryon rest frame is (ai — 4t> m ) cp r 



i e7 



The first term is responsible for the radiation drag due to 
the difference between the wave's material velocity and the 
radiation "velocity" . The second term is due to the nonlin- 
ear effect that all other waves induce on this particular one 
through the change in the medium's opacity and the shear 
between the material velocity and the effective radiation ve- 
locity of all other waves. The value of the large scale shear 
v m o — v r o is the value at the region in which the small scale 
wave is solved for. The latter contribution is synchronised 
with the wave through the changes in the opacity. It is of 
second order and it should be integrated over time. Clearly, 
the contribution of waves oscillating at frequencies different 
than the particular wave averages to zero, thus, the only 
net contribution arise from waves that do not oscillate over 
the relevant integration time; that is to say, only the shear of 
the large scale waves contributes. If a wave oscillates n-times 
during the period, its net contribution will on average fall as 
n _1 . When integrated numerically in §5, the contribution of 
smaller scales is taken into account automatically and their 
contribution indeed falls with uj. 

The radiation field is not perturbed for optically thin, 
isothermal fluctuations. Consequently, v rl vanishes. If we 
define the large scale shear as: 



Av = v m0 
and write 



'<=70 



V r0 , 



po T 



then the radiation term is 

frad — 



Vl | Avp T e71 
7"e7o ^~e7n 7"e 



'=70 ,e T0 



Vl | Avp 



Po T 



(19) 



(20) 



(21) 



The second term does not vanish if the opacity is perturbed 
by the isothermal waves. By definition, the temperature isn't 
perturbed and T\ = 0. If however the opacity is a function 
of the density, then a / and a new term is added to 
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the Newtonian equations of motion. The opacity is predom- 
inantly Thomson scattering, consequently, r is proportional 
to the number of free electrons. The number of free elec- 
trons varies only during recombination. At this time n e is 
inversely proportional to the root of the density. To be more 
exact, the fraction of ionised protons at equilibrium X eq is 
given by Saha's equation: 



damped wave with a damping time scale of r 



The 



- X eq _ 4V2C(3) 



Xeq 



TV \m e J 



3/2 



exp(S ion /T), 



(22) 



with the baryon to photon ratio r\ being proportional to the 
density and Eton the ionisation energy of hydrogen. Through 
differentiation we find: 



d In X eq 
d\np 



1- X et 

2- Xe 



(23) 



That is, before recombination commences, X e? is approxi- 
mately unity and the parameter a vanishes. Afterwards it 
will grow in size and reach the asymptotic value of —1/2. 
Note that for the above equation to be valid, we need 
to assume that the time scale for ionisation equilibrium 

(e + p < > H + 7) is shorter than the dynamic time scale. 

This implies that among other assumptions, the equations 
are valid only before freeze-in which takes place not too long 
after the advent of recombination. We have also neglected 
Helium ionisation and its influence on the density depen- 
dence of the opacity. However, the low abundance of Helium 
results with a much smaller effective a, and together with 
the larger speed of sound, they both reduce the effect during 
the period of Helium recombination. 

Define 8 = pi/po and spatially expand 8, vi and (f>i 
in Fourier integrals proportional to exp (— ik • r/a(t)). The 
result is: 

= 8 k • v fc 



= 



av fc + dv fe — ikv 3 8 k - ik(f> k + a 

4-nGpo 2 c 
<Pk H — a °k - 



Vi 



Av 
aa o k 



1 e~/ 



(24) 



The new phenomena discussed here arise from the introduc- 
tion of the last term in the second equation. The term is 
proportional to the density perturbations (through opacity 
changes) and to the integrated large scale shear Avo which 
varies in space. Even with the additional radiation term, the 
irrotational flow is still decoupled from the rotational part. 
On the other hand, the rotational part is affected from the 
irrotational one. This interesting effect can introduce a seed 
magnetic field from any isothermal fluctuation present, as is 
described in Shaviv & Levin (1998). 

By taking the irrotational part of the second equation 
(through the dot product with k), and using the first and 
last equations, we find that 8 k obeys the equation: 

k+(- + —)s k +( k\l - 4nG P0 - ia^±) S k = 

(25) 

This equation (and its derivation) should be compared 
with the standard derivation of the equations describing 
isothermal waves when all nonlinear interactions are ne- 
glected (e.g., Coles and Lucchin, 1995). The solution of the 
linear equation consists of two modes. One describes a highly 



other describes a frozen mode. Although the latter is for- 
mally unstable, the large radiation drag effectively reduces 
the growth rate to a value much smaller than the universe 
expansion rate. 



4 APPROXIMATE SOLUTION 

The general solution to eq. ^ depends on the exact time 
dependence of a, a, po, Avo & r e7 . We are however, inter- 
ested in the evolution over the short period of decoupling, 
a period that is much smaller than the age of the Universe 
at the time. Moreover, we are particularly interested in cir- 
cumstances were the growth rates are large, in fact, larger 
than the rate in which the various parameters change during 
recombination. This will allow us to assume that the above 
variables are quasi-static. 

The system described by eq. has five different time 
scales of which some depend on the wave number. The first 
two are the expansion age of the universe and the related 
Jeans timescale. For a flat matter dominated universe, the 
first is: 

(26) 

while the second is: 



(27) 



The last equivalence assumed a flat matter dominated uni- 
verse. The third timescale is the isothermal oscillation time: 



\k\v s 



(28) 



The fourth is the electron-photon relaxation time r e7 that 
describes how fast baryon and radiation exchange energy. 
The newly introduced timescale is the time associated with 
the radiation inhomogeneity: 



\J aAvo ■ k/r e . 



(29) 



We will assume here that the electron-photon relax- 
ation time is much smaller than the universe's expansion 
time scale. The reason being that when the timescale sur- 
passes the expansion scale, the solution is in any case invalid 
because the timescale for establishing ionisation equilibrium 
becomes very long. The a/a term can therefore be neglected. 
The time scales themselves depend on the cosmic time as 
well, however, since we are interested in timescales shorter 
than both the expansion of the universe and the timescale 
for recombination of hydrogen, we can for simplicity assume 
them to be constant or quasi-stationary. In such a case, the 
equation's solution is exponential with the roots: 



2r e . 



1 



(30) 



The waves' growth rate is given by the real part of r, namely 
3£(r). It will be governed by the largest rate (or smallest time 
scale) in the system. The oscillatory behaviour is given by 
the Imaginary part - 3(r). For an optically thick universe 
with t c1 < t u , we find two general cases. The two cases 
and their various regions of k in which the rates behave 
differently, are summarised in table |l| and exemplified in 
figure ^ where a few sample graphs of 5R(r(fc)) are shown. 
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Case A: At>o > v s /a 



Boundary between regions 



Growth Rate: R(r) 



T 2 y 

rad 



T rad — T o 



a Avq Tj 



k = 



k k — oc 


!R(r) f 


a t^/t 2 




fc_oo <C k <C fco 


S(r)f= 


- T e-y/ T r ad 


T e7 («Av ■ k) 2 


fc <C fc <C fcoo 


R(r)c 


d ± 1 


I / a Avq -k 
- ±Y 2 r c7 



Boundary between regions 



Case B: Avq < Us/cr 



Growth Rate: 3?(r) 



"7~o — 


fci 


= 1/{t u v s ) — 


Tcy — To 


fc 2 


= l/(r el v s ) — 



k < fci 



fci <C A; < fc 2 



fc 2 < k 



»(r) « t £7 /t 2 



5ft(r) -r e7 /r 2 = --r e7 u 2 fc 2 



Ji(r) Ri -l/2r e . 



Table 1. The growth rate of the small scale waves for various wavelengths under the two relevant cases of At>o > v s /a and Avq < v s /a, 
with T e7 < t u ( i.e., as long as the universe is optically thick). Note that only if SR(r) > t\7 is the solution appropriate. Otherwise, 
instead of an exponential solution, the system will evolve through a power law behaviour with a typical growth scale of r ~ 5R(r). 



When analysing eq. one finds that it admits two 
general forms of behaviour depending on whether Avo is 
smaller or larger than v s /a. The main difference between 
the two cases is that large k waves in the former case are 
damped while in the latter they are unstable and amplify. 
More specifically, one has: 
Case A: Awo < v s /a 

When the shear velocity between the radiation and baryon 
fluid is smaller than the critical value, the shear cannot over- 
come the radiation drag. Even though the quantitative anal- 
ysis may change, the general behaviour is the one found 
when the shear velocity is altogether neglected. That is to 
say, long wavelengths for which t d > r u are frozen in the 
plasma. Although they do have a positive growth rate, it 
is small when compared with the growth rate of the hori- 
zon. Waves with smaller wavelengths for which t < r u are 
damped, and when r e ^ ~ r , the damping rate is saturated 
at its asymptotic value of l/2r e7 . 
Case B: Avo > v s /a 

The behaviour of the system is radically different when the 
large scale shear velocity is larger than the critical value. 
At large wavelengths, the waves are still frozen in and the 
behaviour is the same as in the previous case. On the other 
hand, at shorter wavelengths the waves' growth rate is dom- 
inated by T ra d- Instead of being damped or frozen-in, the 
waves can in this case be unstable. The growth rate of 
waves propagating in the direction of the shear increases 
with k and it saturates at small wavelengths when acoustic 
oscillations become important. The maximum or asymptotic 
growth rate is obtained for values of k larger than: 

fcoo= Av £ ^ i (31) 

i.e., for scales smaller than the Jeans scale by a factor of 
order r„/r e7 . The typical values of r u /r e7 when the rate is 
large are of order 10, implying that the typical scales would 



be 10 times smaller than the Jeans scale. The maximum 
growth rate is then given by: 



R(r) w r /(2r r 2 ad ) 



a Avo 
2 ^7 



n 1 



Te.'y 



(32) 



It is now evident why isothermal waves are crucial. The 
speed of sound v„ in isothermal waves is much smaller than 
the speed of light, therefore, the typical shears that cause 
nonlinear interactions are of order 0(v s /c) and not 0(1)! 
An appreciable nonlinear effect can take place already at a 
dimensionless amplitude of Sp/p ~ 10~ 5 or 10 -4 ! 

We are interested in the integrated effect of the growth 
rate in the time interval between a given to and a given later 
time t\. Hence, we define the following growth parameter: 



G(k,n) 



K(r(fc,n))dt. 



(33) 



If the rate given by eq. [32| is larger than all other inverse 
timescales of the system on which the value of the rate 
can change (i.e., both the expansion and the recombination 
timescales) then the growth factor of the isothermal waves 
is approximately given by: 



c5 tl (fc,n) 



exp(G(fc, n)). 



(34) 



The recombination time At = ti—to is the shortest timescale 
on which the growth rate changes. As a consequence, the 
condition for the approximation to hold can be written as 
5ft(r-) 3> At" 1 , and it can be translated through eq. |33| to 
G> 1. 

We proceed to investigate the rate of growth in two 
ways: 

• The nonlinear growth rate equation (eq. ^) is averaged 
over the distribution of perturbation amplitudes and direc- 
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Figure 2. The growth rate as a function of the wave number depicted for three examples. In all three examples shown v s = 10 c 
and a = 1/2. In example 1: Ano = 3 X 10~ 4 and r e7 = 0.01t u , in example 2: Ai>o = 1 X 10 — 4 and r e7 = 0.1t u , while in example 3: 
Auo = 1 X 10 — 5 and T e7 = 1/30t u . The first two examples correspond to case A of table 1 while the third example for which Svq < v B /a 
corresponds to case B in the table. Only in the former two cases do we find that the growth rate is increased tremendously for large 
values of k. The curve of example 1 has the approximate behaviour corresponding to the four different regions of case A marked with 
small dotted lines in addition to the exact value of the root. In the four different regions, we have k oc const, k 2 , fc 1 / 2 , const. They do not 
exist in example 3, and while they do exist in the second example, the behaviour is not as clear as in the first one. 



tion of motion in order to get the average growth of the 
amplitude of the waves. 

• We explore systems that reach nonlinear amplitudes by 
examining the volume fraction that does develop nonlinear- 
ities and the possible effects triggered by these nonlinear 
regions. 



4.1 Average Growth 

We now proceed to estimate the average growth of amplitude 
of small scale waves for which k <; koo . Although the physical 
scale is small enough not to be directly detectable in the 
near future, these modes have the fastest growth rate and 
are therefore most interesting. 

The av erag ing procedure of the amplification described 
by equation ^ should be divided into two separate averages 
over two distributions. First, the expression should be aver- 
aged over the various possible propagation directions since 
unlike the nonlinear case, the dispersion relation is now a 
function of the direction of k (with respect to the direction 
of the large scale shear) and not only its absolute value. Sec- 
ond, the growth function is found for a given local value of 
the shear, however, it is certainly not linear in its amplitude 
and one should therefore average over the distribution of 
large scale shears. 



We first average expression |34] over the An possible di- 
rections of n (the direction of k) while keeping in mind that 
G is linearly proportional n ■ no. We get: 



St 



Mn) 

5t 



= / dnexp(Goo(fc)(fi • no)) = 



(35) 
sinh(Goo(fc)) 

Goo(fc) 



with Goo(fc) = G(k, n = no), i.e., its maximum value. 

Next, we proceed to average over the distribution of 
shears Avo. This distribution will clearly depend on the 
form of the density distribution. The natural choice for the 
latter is a Gaussian distribution. Although it does not nec- 
essarily have to be the case, it is the theoretically most rea- 
sonable distribution expected for the density fluctuations^. 
The distribution of the absolute value of the shears is simi- 
lar in form to the distribution of the absolute velocities cal- 
culated from density perturbations in large scale structure 
(e.g. Coles and Lucchin §18.3), that is, if the distribution for 
the density fluctuations is Gaussian then it will be Guassian 



* Although it is a natural consequence of inflationary scenarios, 
it is markably different for models in which cosmic defects provide 
the origin of the fluctuations (e.g. Turok 1996). 
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for the components of the shear vector in each axis, and 
the distribution for the absolute value will be given by a 
Maxwellian distribution. If the root mean square of Goo is 
Grms, the distribution itself is: 



P(G 



54 Gl 



■ exp 



/ Goo y 

\ Grms J 



7T Gl 

The average perturbation growth factor is therefore 



(36) 



St, 



/ all 




exp 



G r 



(i 



exp 



A(Gr 



f Goo y 

V Grms / 



sinh(G c 

Gog 



(37) 



-dGo 



The function _4(G rms ) describes the amplification fac- 
tor of the fluctuations. Clearly, for Grms much less than or 
of order unity, the effect is small or negligible. Moreover, the 
numerical values obtained, such as a 10% increase for G = 1, 
are only approximate since the required condition G>1 for 
the validity of eq. is not fulfilled. Nonetheless, already at 
Grms = 5, small scale isothermal fluctuations will amplify 
by almost 2 orders of magnitude. They will grow by 5 orders 
of magnitude for G rms = 8 and for G rms = 12 the growth 
will be a staggering 10 orders of magnitude. Note that if the 
growth is saturated when the amplitude reaches nonlinearity 
then the average growth is meaningful only if the amplified 
amplitude is less than unity. Moreover, unlike the case when 
Grms ~ 1, the solution here is much more accurate because 
regions that contribute the most to the growth are those in 
the high Goo tail, where the approximate solution is valid. 
For example, the average growth when G rms = 8 is of 5 
orders of magnitude and the main contribution originates 
from Goo ~ ln(10 5 ) w 11.5 > 1. 

Also worth mentioning is that in the averaging proce- 
dure over the fluctuation distribution, we are essentially av- 
eraging an exponential function ~ exp(Goo); thus, if the 
distribution for the amplitudes does not fall fast enough for 
large values of Goo (which is proportional to Sp/p), the aver- 
aging integral diverges. Such a divergence occurs if the high 
5p/p tail is for example a power law. The implication is that 
there will always be a finite volume (i.e., volumes that are 
not exponentially small) where nonlinearity can be reached 
if the distribution has a wide tail. A wide tail can in this 
case be any distribution with a power law tail or even a dis- 
tribution with an exponential tail, provided the tail of the 
distribution is in the form oc exp (—Goo /Go) with Go > 1. 



4.2 The nonlinear regions 

The large exponential growth associated with small regions 
where the shear is large, is of high importance. In these re- 
gions that are small for a modest value for G rms , the growth 
may be large enough as to amplify waves out of the linear 
regime and let the nonlinear evolution of these wave com- 
mence. Consequently, it is interesting to estimate what frac- 
tion of space will actually be filled with nonlinear waves by 
the end of recombination. 

The largest growth factor in a given region is given 



by ~ exp(Goo). Thus, if the typical small scale fluctuation 
amplitude before amplification is initially (5o, then only if 
5oexp(Goo) iJ 1 will there be waves that grow out of lin- 
earity, i.e., only if Goo ii —In So- If the density fluctua- 
tions are Gaussian, the distribution of growth factors Goo 
is Maxwellian and given by eq. |3(| The fraction of space in 
which there are waves that reach nonlinearity is then given 
by 



F{8o, G r 



54 Gl 



In S 



7T G'rri 



■ exp 



3 / Goo y 

2 V Grms ' 



dG 



= erfc 



\ A KJrms / V \ r m s 



In So 

Grms 



(38) 



Here erfc denotes the complementary error function. The 
function T is plotted in figure ^| As an example, for G rms ~ 
5, So should be roughly 10 -4 to have 1% of the volume reach 
nonlinearity. 

If we could observationally detect all the regions in the 
background microwave sky that reach nonlinearity (see e.g. 
§5) then the absolute number of regions showing nonlinearity 
will be given by the fraction that does develop it multiplied 
by the total number of statistically independent regions. The 
latter number is roughly 4ty£ 2 with £ the multipole order cor- 
responding to the typical shear producing wavelengths. By 
taking a co- moving k ~ 0.1A/pc _1 as the typical correlation 
length for the shear producing zones (see e.g. §5 or fig. 
we get: 



£ w (6000 h~ 1 Mpc)k ~ 1000, 



(39) 



namely, to get a statistically significant number of regions, 
one must have: 



F(5q, Grms) » 



1 



4n£ 2 



10" 



(40) 



Evidently, for our example of G ~ 5, any initial isothermal 
perturbations that has So <; 10~ 6 is theoretically detectable 
and will produce at least 100 different nonlinear regions. If 
the primordial spectrum of fluctuations is isothermal and 
flat, and the expectation for So is ~ 10 -4 then for any 
Grms > 3 one should in future be able to detect the effect. 
For G <; 10 more than 50% of the universe reaches nonlin- 
earity at 2 ~ 1000. In the case where the isothermal wave 
amplitude has the smallest plausible value, namely when 
the primordial spectrum is flat and adiabatic and the Press- 
Vishniac (1979) effect produces isothermal waves with an 



amplitude do ~ 10 
any Grms ^ 7.5. 



10 , the effect will be detectable for 



5 ESTIMATING THE EFFECT IN ACTUAL 
COSMOLOGICAL SCENARIOS 

The largest attainable growth parameter is roughly: 
a At; At 

2 Us Tfi'V 



G 



(41) 



where At is the approximate duration of decoupling. v s is 
the isothermal sound speed which at recombination is about 
4.8 orders of magnitude less than the speed of light. There- 
fore, when recombination commences and r e7 ~ 10~ 2 r u (for 
an approximate duration of AT w 10 _1 r u ), the shear veloc- 
ity should be about 3 x 10 c to get a 2 order of magnitude 
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Figure 3. The fraction of regions that reach nonlinearity — JF(<5o, G rms ) as a function of Logio(l/So)/G rrns . Any fraction less than 
~ 10~ 6 will be theoretically impossible to detect and correlate with gradients in the CMBR at a 3cr level, as there are only about 10 7 
statistically independent regions for which the shear is substantially different and uncorrelated. 



increase and 6 x 10 _5 c to get as much as a 10 orders of 
magnitude increase. As recombination continues, the shear 
increases while l/r e7 and the optical depth decrease roughly 
by the same factor, such that the rate of growth is effec- 
tively unchanged during this period, even though the opac- 
ity changes by more than two orders of magnitude. 

The interesting question that will ultimately determine 
the fate of the small scale fluctuations is therefore a quanti- 
tative one - what is the exact value (i.e. the exact r.m.s.) of 
the shear exerted by large scales and what is the integrated 
rate equal to? 

To answer the quantitative question, a numerical simu- 
lation for the behaviour of the large scale modes was carried 
out. The numerical code is a modification of the COSMICS 
code developed by Bertschinger and Ma. The code was used 
to integrate the linearised equations of general relativity, 
matter and radiation in the synchronous guage, and sub- 
sequently normalize the power spectra to the COBE mea- 
sured quadrupole moment. A description of the equations 
solved and the numerical algorithm used can be found in 
Bertschinger (1995) and Ma & Bertschinger (1997). With 
the modification, the maximum growth rate given by eq. ^ 
is calculated and integrated to give the growth parameter. 

The evolution of each co-moving large scale wave vector 
k is followed over the course of the particular scenario's cos- 
mological history. For each k integrated, the code calculates 
the evolution of the wave amplitude of the various species 
in the primordial soup, including among other, the angular 
moments of the photon distribution Ae(k, t). These are to be 



precise, the Fourier and Legendre decomposition of the pho- 
ton brightness temperature perturbation A(x, n, i) = AT /T 
of the photons traveling in the direction of n; more specifi- 
cally, they are defined through: 

/oo 
d 3 ke^ x ^2(-l) e (2£+i)A e (k,t)P e {k-n). (42) 

1=0 

The program also integrates the appropriate unnormalized 
growth factor: 

,, . f a Avn dt 
S»(M)= / o— — . ( 43 ) 
Jo Vs Te ~y 

where Avq is the unnormalized shear. 

To normalize both Ag(k,t) and goo{k,t), one should 
take into account the original, primordial power spectrum 
P(k) = Ak 4 ~ n . When n = 1 the power index corresponds 
to a flat Harrison-Zel'dovich power spectrum for which the 
total power per decade is constant. We normalize the spec- 
trum to fit the measured COBE quadrupole moment. This 
is achieved through the separation of the angular correlation 
function C(9) into its moments: 

C(6) = J2 ^p-C e P e (cose), (44) 

e 

with Pt the Fth Legendre polynomial, and writing the mo- 
ments as a function of the normalized photon power distri- 
bution functions: 
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Run 


Model 


n b 


n c 




n 


Grms 


1 


std. CDM 


.05 


.95 




1 


5.1 


2 


std. CDM 


.02 


.98 




1 


4.7 


3 


tilted. CDM 


.05 


.95 




.9 


4.2 


4 


Mixed HCDM 


.05 


.75 


.2 


1 


0.9 


5 


A-CDM 


.05 


.55 


.4 


1 


2.8 



Table 2. Summary of the numerical results. All models are 
flat (curvature free) universes with a Hubble constant of Hq = 
66fcm sec -1 Mpc~ 1 and with adiabatic initial perturbations. 



C £ = 4tt / d 3 kP(k)A 2 e (k). 



(45) 



The quadrupole moment is related to the second moment 
Ci of the angular correlation function by the definition: 



Qr 



T 



1/2 



(46) 



through which the constant A of the power spectrum can 
be normalised. The value used for the quadrupole moment 
is Qrms-ps = 18juisT Q The r.m.s. of Goo can then be re- 
lated to the unnormalized growth rate (eq. ^3|) through the 
expression]^: 



(47) 



A few cosmological scenarios were checked, all of which 
correspond to flat universes. Typical results are summarised 
in figures ^-|^. In figure ^ we find the various physical pa- 
rameters that affect the large scale shear. Figures ^ and ^ 
depict the unnormalized growth rate and its integrated value 
(the growth parameter) as a function of time for several co- 
moving wavenumbers in a standard CDM model (run 1). 
The last figure describes the integrated growth parameter 
as a function of co-moving wavenumber for two scenarios - 
a standard CDM and a mixed Hot-Cold dark matter sce- 
nario. We find that most of the contribution comes from 
a co-moving k ~ 0.1 Afpc -1 but a significant contribution 
arises from scales that are smaller or larger by as much as 
an order of magnitude (or even more in the case of CDM) . 

The scenarios simulated are summarised in table |^. We 
find that the typical r.m.s. of Goo varies from slightly less 
than unity to more than 5, i.e., the average growth of the 
small scale isothermal waves in flat cosmologies can range 
from 10% to 2 orders of magnitude, depending on the par- 
ticular scenario. 



t The four year 1 - a COBE DMR results are Q rms -PS = 18 ± 
1.6 fiK if n is assumed to be 1 (a flat spectrum). If it is not 
constrained, then Q r ms — PS = 15.3^11 ^ K ( Benn et et al. 1996). 
The final value of G r ms is proportional to the value of Qrms— PS 
used. 

t The coefficient before the integral (unity) is chosen to be con- 
sistent with P(k) defined to be the total power spectrum, thus 
consistent with eq. W5 as well. 



6 CAN EARLY SMALL SCALE STRUCTURE 
FORMATION BE MEASURED? 

We have found that small scale perturbations can under 
certain conditions be enormously amplified. The strong en- 
hancement of particular waves leads to the following ques- 
tions: First, if the required conditions do prevail and the 
amplification factor is more than enough to drive the per- 
turbation out of the linear regime, at what point will the 
growth process stop? Second, if this process did in fact take 
place early in our universe, could it have left a mark on the 
CMBR that is detectable today? And third, could the pro- 
cess have had other implications on structure formation in 
the early universe? 

It is difficult to precisely estimate when the growth sat- 
urates. The simplest estimate to when the growth process 
ceases is as the wave reaches the nonlinear regime, a point 
at which our linear treatment breaks down anyway. Clearly, 
when 5p/p is of order unity, the different speeds of sound at 
different parts of the wave can create shocks that will dissi- 
pate the waves' acoustic energy into thermal energy. When 
8p/p ~ 1, the acoustic energy of isothermal waves is of or- 
der the thermal energy of the gas alone (In adiabatic waves 
the energy is of order the thermal energy of the gas and 
radiation). At the time of recombination, the heat capacity 
is very large because a large amount of energy is needed 
to reionise the Hydrogen. Thus, after recombination, when 
the temperature is of order but less than the recombination 
temperature T & T rec w 0.26eV, injecting an amount of 
order the thermal energy of the plasma will ionise at most 
only ~ T rec /E io „ « 0.26el//13.6eV « 1/50 of the Hydrogen 
atoms. The free electrons released by ionization rescatter 
photons and leave a mark in the CMBR. 

Heating the plasma back to the recombination tempera- 
ture implies that the plasma and the microwave background 
are not in equilibrium anymore. The plasma consequently 
upsets the Planck distribution of the radiation through the 
number conserving Compton scattering process. The pro- 
cess is evident from Kompaneets' equation (1957) and its 
cosmological application by Sunyaev and Zel'dovich (1969). 
The dimensionless parameter describing the importance of 
the effect is the parameter y defined as: 



fc(T e > 



(48) 



with standard notations. n e denotes the number of free elec- 
trons. The analysis of the COBE data revealed no deviation 
from a Planckian spectrum and set a limit of (Fixsen et al. 

1996) :§ 



y £ 1.5 X 10" J (95%CL). 



(49) 



If the fraction of ionised electrons increases to i e ~ 1/50, 
and the temperature T increases to tT with r ~ 1, one can 
show that the resulting y parameter should be of order (e.g. 
Peebles, 1993): 



§ Note that a smaller upper limit of y Si 3 X 10 -6 (Fixsen et al. 
1997) was found when the signal was correlated with the DMR 
fluctuation map, however, this type of limit assumes that y is 
correlated with the temperature fluctuations, in our case however, 
the correlation will be with the temperature gradients. 
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Figure 4. The various physical variables that affect the net shear in the system and its effect on small scales as a function of expansion 
radius of the universe for a sCDM model with Hg = 66 km s~ 1 Mpc~ 1 and = 0.05. The variables are a = d log X e /d log p which 
directly affects the strength of the force induced on small scales, X e , the number of free electrons to which the total force is proportional, 
t/rej gives the mean free path of an electron in the photon "fluid" to horizon scale ratio, while X e /X e fi na i is the ionisation to residual 
ionisation ratio. If it is much greater than unity, the reaction e + p < > H + 7 is still in equilibrium and the solution is valid . 



y ~ 2 x 10 _11 (r- 1)(1 + z) 5/2 hfl b n- 1/2 L e . (50) 
with z ~ 1000 at recombination. Thus, 

Umax ~4X 10~ 6 /lo. 75^60. 05^ _1/2 - (51) 

If through the appropriate nonlinear treatment it would be 
found that the kinetic energy of the saturated waves is larger 
than the thermal energy then the energy, disspated is larger 
and y m ax should be increased correspondingly. If on the 
other hand the dissipation is found to take place on smaller 
time scales than the universe's age at the time, r — 1 is 
smaller and it reduces y m ax- The actual values cannot be 
found from our linear analysis. 

The aforementioned value for the y parameter is how- 
ever only in regions (and directions in the sky) where the 
gradient at recombination was large enough to have isother- 
mal waves reach saturation. The average y will be smaller 
and given by: 

„ -y rf 1 ^^) (52) 

yavr — ymax-J I -, . , I ■ \va j 

\ log G y 

Unless the fraction T of the sky that reaches nonlinear- 
ity is of order unity, the y avr is too small to be detectable by 
present means. Nevertheless, investigations that correlate y 
with the gradients in future sub-degree maps will be able to 
detect smaller values for y. 

If the observations rule out the existence of such struc- 
ture then the theory imposes strict limits on the cosmo- 



logical parameters, including the fact that the fluctuation 
distribution function does not have a power law tail. How- 
ever, it can predict black-body deformations correlated with 
the CMBR gradients (i.e., a specific correlation). If indeed 
this is the case, it will be proven as the origin for small scale 
structure, and it could be used for setting strict limits on 
cosmological parameters. 



7 EXOTIC POSSIBILITIES 

More exotic possibilities for the fate of the nonlinear re- 
gions rely on the fact that the sub- Jeans sized objects formed 
through the radiative instability during recombination can 
later on gravitationally collapse. The typical wave numbers 
that experience the large amplification are those with a 
wavelength larger than the Jeans scale by at least an or- 
der of magnitude. If they exist in the linear regime they are 
unstable and collapse when z ~ 10, however, these wave 
lengthes are already nonlinear at z ~ 1000 and may eas- 
ily collide and collapse long before their linear cousins will. 
Compact objects of order of the Jeans scale can therefore 
form and collapse at very high redshifts 10 Si z & 1000. 
What will the fate of these objects be? They might collapse 
into massive Black holes, they can also burn and release a 
fewxl0 _J of their energy while doing so (Bond et al. 1984). 
These nonlinearities may perhaps form globular cluster like 
clusters or maybe even seeds for AGN engines. 
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Figure 5. The (unnormalized) evolution of the instability growth rate g oc (k,t) as a function of the universe expansion scale for five 
different wave numbers; namely, the contribution to the large scale shear from various wavelengths during the time of recombination. 
For sub-Horizon scales which are not too small, waves with larger k oscillate more times during recombination and cancel out some of 
their contribution. However, when k is increased further, the optical depth corresponding to one wavelength falls below unity and the 
radiation fluid and baryon fluid stop oscillating synchronously. As the baryon fluid loses its radiation support, its oscillations decreases, 
as can be seen in the figure. As a consequence, the contribution from large wavenumbers does not fall as fast as 1/fc, as is evident in fig. 
8. The scenario is a standard CDM with Hq = 66 km s~ 1 Afpc -1 , f2(, = 0.05 and a flat adiabatic primordial spectrum (run 1). 



Another possibility is that the nonlinear objects form 
structure on larger scales through the large energy they re- 
lease and a chain reaction of explosions that they ignite (Os- 
triker & Cowie, 1981), they can then prompt the collapse of 
galaxies in the expanding shells and form voids. Although 
the peculiar velocity and the y-distortion prediced for the 
simplest of these models are both too high (Levin, Freese & 
Spergel, 1992), the inclusion of continuous or multi-cycled 
detonation waves in the explosion model can lead to large 
voids with a diameter larger than 10 Mpc, a peculiar ve- 
locity of roughly 100 km sec -1 and a ^-distortion less than 
the limit set by COBE (Miranda & Opher, 1997). If these 
explosions can indeed be ignited, only a very small fraction 
of the universe needs to reach the nonlinear regime by the 
end of recombination in order to have dramatic effects on 
the formation of large scale structure. 



8 FEEDBACK ON LARGER SCALES 

The next question to address is whether the effect can also 
affect the large scales. Since energy is conserved, large scales 
will be affected if the energy transferred to the smaller scales 
is comparable to the original energy in the large scales. The 
energy per unit mass in the large scale waves is roughly 
given by (Sp a d/p) 2 v 2 a - j /2 while the energy per unit mass in 



the small scale isothermal waves is roughly (5pi so / p) 2 v 2 i - ) /2. 
Inspection of eq. |l^ reveals that when the small isothermal 
waves have a dimensionless amplitude of order unity over 
a large fraction of the universe and G rm s ~ 1, the energy 
content in both types of waves is comparable. A significant 
fraction of the energy in the large scales can therefore be 
transferred to the small scales. 

Apparently, an exact treatment is beyond the scope of 
the linear analysis developed here and the unknown nonlin- 
ear physics can at most be phenomenalised using few pa- 
rameters and a simple demonstrative model. 

Using the terminology and approximations of §2, one 
can express the original energy per unit mass in the large 
scales as A 2 G 2 ms x v? i \/2. Under the approximation of §2, 
we have Ai — 4, however, for optically thin large scale waves 
the energy in the waves is different and Ai can change. The 
energy in the isothermal waves is A 2 x v 2 i - j /2 and the energy 
growth rate is 2G r msA 2 /At x v?^/2. 

Next, we assume the energy damping rate of the small 
scale waves to be of the form: 

u <- damping 

At is the duration of recombination. Ad and e are parameters 
that represent the unknown physics. Both parameters can 
depend on the wavelength of the unstable isothermal waves. 
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Figure 6. The same as the previous figure except that plotted are the time integrated growth rates. 
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Figure 7. The (unnormalizcd) integrated growth rate as a function the co-moving wavenumber k for two cases. One is the standard 
CDM (run 1) while the other is a mixed Hot-Cold dark matter (run 4), with parameters given in table 2. 
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Ad should be of order unity for waves with a period of ~ 
At, but otherwise unknown. The power e with which the 
damping depends on the amplitude of the isothermal wave 
is also unknown but should satisfy e > 2 if it is characterised 
with only one power index and if it is to be more important 
on the nonlinear scales than on the linear ones. 

If the value of G r ms before the instability commences is 
large enough then a large fraction of the volume leaves the 
nonlinear regime. The critical value for this to take place is 



Grms <; — hx(5piso/p)' 



(54) 



The simplest simplifying assumption is that the isothermal 
waves will saturate at a level where the energy growth rate 
is equal to the dissipation rate. The small scale amplitude 
A„ will then be given by 



main conclusion is that for a large range of an initial large- 
scale amplitudes, it is possible to have almost all the energy 
dissipated into smaller scales and subsequently have it dis- 
sipated through nonlinear processes. The large scales then 
reach their natural lower limit - the value for which the large 
scale shear switches off the instability. 

The implications are very interesting. The measured 
fluctuations in the microwave background radiation reflect 
the large scale structure after the fluctuations were repro- 
cessed by the instability during recombination. In other 
words, the value oberved does not necessarily reflect the 
pre-recombination value of the fluctuations. It is possible 
to have an initially large large-scale amplitude, induce non- 
linear small scale structure and eventually reduce the large 
scale amplitude. 



A 3 = 



(Grms \ 
~Aa~J 



(55) 



Compare now the energy in the large scales with the total 
energy that can be dissipated. When the latter energy is 
larger, all the energy in the large scales can be dissipated. 
The condition for it to take place is: 



G 

Grms ^ 2 



< 2 t - 2 A 2/{4, ~ t) A 2,(A ~ e) 



A 

rms rO ^ d 

2^2/(4-e)^2/(4-e) 



for e < 2 or e > 4 
for 2 < e < 4 



(56) 



When 2 < e < 4, an initial value of G rm s that satisfies 
the above condition will imply that the amplitude of the 
large scale fluctuation can be reduced until G rm s does not 
comply with the above condition. For other values of e, the 
initial value of Grms must be smaller than a critical value in 
order to satisfy the above condition, in which case, the large 
scale amplitude will apparently be reduced arbitrarily. In 
fact, Grms continues to decrease until the shear drops below 
the critical value separating between the two cases, at which 
moment the instability turns off. Since r e7 ~ 1 at the end of 
recombination (see e.g. figure^), the switching off will take 
place at G rms ~ 1. To summarise, two possibilities for the 
damping of large scale waves seem to exist depending on the 
power e relating the dependance of the damping rate with 
the amplitude of the isothermal wave. In the first case, a ny 
Grms larger than roughly the critical value given by eq. |56J 
will be reduced to a value of order of critical value provided 
that it satisfied eq. as well. In the second case, any Grms 
which satisfies both eq. ^6] and eq. |f34|, will be reduced to 

Grms **** 1- 

As an elucidation, one can numerically solve a simplified 
toy model for the process. The transfer of energy from the 
large scales into the small scales and then the dissipation 
by small scale nonlinearities can be simplified and described 
with the two differential equations: 



d 



— (j4f Grms(i) 2 ) — — 



+ - 



2>Grms (t) 

AT 

2GV 



As{tf 



(57) 



Figs. |8ja & |8p describe the numerical results for the integra- 
tion of the two differential equations when e = 3, Ad = 1, 
Ax = 4, for initial fluctuations of Spiso/p = 10 -5 and differ- 
ent initial values of G rm s- By comparing to the numerical 
results to eq. it is clear that the analytical estimate is 
an over simplification of the differential equations, which 
too are a simplification of the real physics. Nevertheless, the 



9 DISCUSSION 

We have analysed the nonlinear growth of small scale 
isothermal fluctuation through the induced effect that large- 
scale perturbations have. The equations are solved through 
the separation of the large scale inhomogeneities into smaller 
regions where they can be considered uniform, and the shear 
Avo exerted by the large scales is constant in space. This can 
be performed because the waves that interest us are waves 
with very small wavelengths and a very small propagation 
speed, such that over the integration time, the waves prac- 
tically remain in the same region of space. 

Small scale isothermal waves are found to be unsta- 
ble if the large scale shear velocity between the "photon 
fluid" and the baryon fluid is larger than roughly twice the 
isothermal speed of sound at the time of recombination. A 
period when the logarithmic derivative of the opacity with 
respect to the density does not vanish but is —1/2 instead. 
When the shear is smaller than the critical value, the solu- 
tion found is qualitatively the same as the solution for the 
completely decoupled equations. The shear in this case may 
induce some quantitative corrections but the behaviour of 
the system does not change. The behaviour of the system 
is however radically different if the shear is larger than a 
critical value. For wavelength shorter than the roughly the 
Jeans scale, the growth rate increases with k and it saturates 
when r e7 = t , i.e., at a scale which is a few times smaller 
than the Jeans scale. The maximal growth rate attained is: 



(58) 



_ a Avo 1 

°° 2 V s T e j' 

If Too is much larger than all other rates in the system, it 
can be easily integrated over the duration of decoupling At. 
The growth factor obtained is then: 



5i , . ( a Avq 1 
— w exp(Gcx)) = exp ( 

do V \ 2 V s Te7 



At 



(59) 



with the brackets denoting a temporal average. The isother- 
mal waves are unstable because their sound speed is very 
small. Consequently, the acoustic energy in a wave with 
a given amplitude is very small and the amount of work 
needed to increase its amplitude isn't large, in fact, it is so 
small that nonlinear effects start to take place already at 
5p/p~ 10" 5 to 1(T 4 ! 

Even more surprising is the fact the the typical shears 
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Figure 8. The two panels represent the results of the numerical integration of the two approximate differential equations describing 
the evolution of the energy in the large scales (A^G^ ms X u^. /2) and in the small scales (A 2 S X /2). The initial small scale amplitude 

before the instability takes place is A s (t = 0) = 10 . Also, e = 3, = 1 and A[ = 4. The evolution of A s (t) and G rms (t) is depicted in 
panel a for various initial values of Grms- Note that only part of the recombination period {0, At} is shown. The same cases are drawn 
with the same line type in both parts. Panel b depicts the final values at the end of recombination - A s (t = At) and G rm s(t = At). 
For the switching between cases A & B of table 1, it was assumed that the critical shear corresponds to G rms = 1. It can however be 
somewhat larger or smaller. 
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found in typical cosmological scenarios is exactly at the verge 
of having nonlinear effects take place! If the amplitude of the 
large scale fluctuations would have been an order of mag- 
nitude larger, then any isothermal perturbations, however 
small, could have been amplified right into the nonlinear 
regime. On the other hand, if the typical amplitude would 
have been a few times smaller, the effect would have been 
completely meaningless. 

For typical cosmological scenarios normalized to COBE, 
the growth parameter is roughly G rm s ~ 1 — 5, e.g., for an 
« = 1, n b = 0.05, h = .66 CDM model one finds G rm3 w 5, 
if the spectral index is tilted to n = 0.9, one finds G rms ~ 4, 
if a hot component is added at a 20% level to the untilted 
model, it falls to G rm s ~ 1. 

The fate of the regions that do reach nonlinearity is still 
an open question. Do these regions collapse and form black 
holes? Do they form an early generation of massive stars? Do 
they release enough energy to the environment and affect it, 
or perhaps, even form structure on much larger scales? Al- 
though some speculations of what might occur do exist in the 
context of early small scale structure formation, it is not en- 
tirely clear. Moreover, nonlinear hydrodynamic simulations 
is probably unavoidable if we are to really solve the problem. 
The only relatively certain consequence is that the dissipa- 
tion taking place after recombination heats the matter and 
it is likely to ionise it, raise its temperature to the ionisa- 
tion temperature and leave an imprint on the CMBR in the 
form of a deviation from a Planckian spectrum. In scenarios 
where G rms <; 5, isothermal perturbations with an ampli- 
tude of as low as 5pi so /p ~ 10~ 6 will be in fact theoretically 
detectable in the future. An amplitude of 8p iso /p ~ 10~ 4 
will leave 1% of the background sky with a detectable devi- 
ation. This fraction can change considerably if the density 
perturbations are not Gaussian. 

The detection of deviations from a Planckian spectrum 
or the placing of an upper limit for them can result with in- 
teresting implications. Any deviations found will first of all 
directly prove the existence of primordial isothermal fluc- 
tuations. Second, such a detection will place stringent lim- 
its on cosmological parameters, as only those scenarios that 
produce a large enough G rm s are capable of producing non- 
linear structure, and because the amount of nonlinearity is 
extremely sensitive to Grms, it is also sensitive to the ex- 
act cosmological parameters. Third, if the cosmological pa- 
rameters are known with an large accuracy (for example, 
through the fitting of the CMBR spectrum) then the pri- 
mordial isothermal perturbation spectrum at very large k's 
can be estimated. 

Even if no detection of a deviation from a Planckian 
spectrum is found, one can place limits on cosmological pa- 
rameters. Moreover, if in the future it will be found that 
the cosmological parameters actually correspond to a high 
G rm s model, one will be able through the lack of detection 
of a y-parameter to place extremely powerful limits on the 
amplitude of the isothermal component of the primordial 
spectrum. Through the Press & Vishniac effect, limits can 
also be placed on the adiabatic spectrum. 

Another interesting implication is the possibility of 
transferring enough energy from the large to the small scales 
and considerably change the amplitude of the large-scales. 
This possibility relies on whether a significant volume frac- 
tion can be amplified out of the linear regime and on the 



maximum dissipation rate of nonlinear isothermal waves. 
Under favourable circumstances, the large scale amplitude 
can be significantly reduced such that G rms calculated from 
the observed large scale fluctuations would only be of or- 
der unity. The fact that the observed fluctuations in the 
CMBR correspond to values of this order raises a very in- 
teresting question. Why do the values of G rm s corresponding 
to the observed fluctuations happen to fall in a small region 
around unity? Is it because the universe was created with 
a fluctuation spectrum corresponding to this region, or, is 
it because G rms that evolved from the initial cosmological 
parameters was actually larger but it was naturally reduced 
to the oberved value? 

One can summarise the several plausible scenarios ac- 
cording to both the value of G rms predicted by the cosmo- 
logical model and the amplitude 5 = 8pi SO /p w 10 -10 — 10~ 5 
of the sub Jeans scale isothermal fluctuations before the ad- 
vent of recombination, these possibilities are: 

(i) If Grms is of order unity or less, the typical growth 
of isothermal waves is at most a few e-folds. The effect will 
be insignificant as the predicted micro degree size fluctua- 
tions are neither fluctuations of the temperature nor on a 
scale measurable in the near future. Moreover, there are no 
implications at all to structure formation. 

(ii) For a value of G rm s iS — (ln<$)/3 ~ 4 — 8, the effect 
will be measurable as patches in the CMBR with a distorted 
Planckian spectrum. As long as G rms Si — In 8 w 10 — 20, 
only a small fraction of the universe would have reached 
nonlinearity and formed small scale structure by the end 
of recombination. Under certain circumstances however, it 
can lead to large scale structure formation through explosive 
amplification. 

(iii) For a value of G rms such that G rm s ii — In 8 « 
10 — 20, a large fraction of the universe reaches nonlinearity 
by the end of recombination and small scale structure is sub- 
sequently formed in most of the universe soon after recombi- 
nation. The large scale spectrum is damped by transferring 
a significant amount of energy to the small scales, thus, the 
Grms measured from the damped spectrum is smaller than 
the Grms calculated when neglecting this process. In some 
cases, Grms will be reduced to a values of order unity and 
it will mimic values on the boundary between the first and 
second scenarios. Note that it will not affect fluctuations in 
components that decouple from the matter-radiation fluid 
before recombination (e.g. dark matter fluctuation). 

The analysis presented here is the first step in the inves- 
tigation of the radiation-matter interaction instability. More 
accurate integration is needed to improve the actual trans- 
fer function for isothermal waves. The analysis here was re- 
stricted to the evaluation of G rms and it does not include 
the actual rates which are k dependent. This approxima- 
tion overestimates the rate of growth of finite sized fc's. In 
the present paper we have simulated the freezing of the in- 
stability due the freezing of recombination by stopping the 
integration abruptly at a given z. However, the equilibrium 
conditions and therefore the switching off of the effect de- 
pend on the isothermal wavelength as well. By assuming this 
assumption, we have actually underestimated the contribu- 
tion from waves of order the Jeans size. 

Many of the possible implications depend on the quan- 
titative behaviour of of the nonlinearities once they are 



reached. A numerical hydro dynamic study of large ampli- 
tude waves will certainly help us understand of the fate of 
the nonlinear objects and the possible implications they have 
on large scale structure as well. 
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